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Abstract —The probability of large-angle scattering for multi- 
GeV muons in thin lead targets (t/Xo of 0(10“^)) is studied. 
The new estimates presented here are based both on simulation 
programs (GEANT4 libraries) and theoretical calculations. In 
order to validate the results provided by simulation, a comparison 
is drawn with experimental data from the literature. This study is 
particularly relevant when applied to muons originating from 
CC Interactions of CNGS beam neutrinos. In that circumstance 
the process under study represents the main background for 
the Ufi Vr search in the r —> /i channel for the OPERA 
experiment at ENGS. Einally we also investigate, in the CNGS 
context, possible contributions from the muon photo-nuclear 
process which might in principle also produce a large-angle muon 
scattering signature in the detector. 

Index Terms —Large-angle muon scattering. Lead, form factor, 
simulation, OPERA, Vr appearance. 


I. Introduction 


M ultiple Coulomb scattering has been studied exten¬ 
sively in the past in many different conditions and for 
different incident particles both theoretically and experimen¬ 
tally. In this work we are interested in the characterisation of 
large-angle Multiple Coulomb Scattering (LAS) of multi-GeV 
muons scattering over high-Z materials with a thickness being 
a fraction of the radiation length Xq (i.e. thickness 0(1 mm)). 
In the LAS regime we are interested in, the transferred mo¬ 
mentum q is such that the De Broglie wavelength (A ~ h,/q) 
of the probe is much smaller (< 1 fm) than the typical radial 
extension of a high-Z nucleus (about 6-7 fm for Lead). For 
a correct description of the high-angle tails of the scattering 
angle distribution it is thus mandatory to precisely take into 
account the effect of the nuclear charge distribution as well as 
the contribution of scattering off individual protons. 

In Sect. |I^ we review basic theoretical aspects of the problem 
highlighting the associated challenges. In Sect.|^we overview 
the available approaches for simulating multiple Coulomb 
scattering while in Sect. |IV] we describe specifically the 
simulation program developed for this work using GEANT4 
(ID, 0) libraries. In Sect. [V] we compare the developed 
simulation with a theoretical approach described in Sect. 

In Sect. VI we validate the simulation on experimental data 
coming from past experiments designed to study the nuclear 
charge density. Finally, in Sect. ivnl the impact of LAS as 
a background for the search is evaluated. The 

importance of muon photo-nuclear interactions in mimicking 


LAS topologies at CNGS is evaluated in Sect. VIII 
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II. Multiple scattering theories 

Multiple scattering has been an active field of theoretical 
and experimental efforts since very early days of particle 
physics. Goudmit-Saunderson E, a and Lewis a developed 
a theory which allows to calculate the exact distribution of 
the angle in space (9) due to multiple scattering after a 
given path s. The problem was tackled by formulating it in 
terms of a transport equation solved by expanding the angular 
dependence in Legendre polinomials (P/(cos0)). The angular 
distribution is found to be given by: 

°° 01 4- t 

F{0,s)=Y^ e~VP;(cos 6 l) (1) 

1=0 

with the transport mean free paths defined as: 

-1 -2TrsN (1-Pi (cos x))^^^d(cosx) 

A; =e -^-1 . ( 2 ) 


The theory does not assume any particular form of the single 
scattering differential cross section da/dXl and provides exact 
solutions. It can also deal with arbitrary nuclear form factors 
provided that integrals are evaluated numerically. Nevertheless 
for large scattering angles the calculation requires the evalua¬ 
tion of high order Legendre polinomials which are difficult to 
treat. This limitation is also relevant for the present case and 
hence this approach was not considered further. 

Another interesting theoretical study was performed in Q 
by Meyer. The theory by Moliere Q, 0 is there extended to 
take into account a generic nuclear form factor. The general 
distribution of the scattering angle 9 is derived and, for a point¬ 
like nucleus, it reads as: 


F{9,t)9d9 = 2xe-^" 


1 + 


bp + b2X^ + EIL 2 
B 


2i/-| 


dx 


(3) 

X being equal to 9 apart for numerical constants which depend 
on the parameters of the problem (incoming momentum, 
particle, target thickness, density etc.) which are made explicit 
in Appendix A. The other quantities are either constants 
&o = 0.423, 62 = —1.423, & 21 / = simple functions 

of the system parameters as for B (see Appendix A for more 
details). The component having the form x^~^ results from 
multiple consecutive scattering in the target thickness while 
the tail at higher x arises from the occurrence of single high 
angle scattering. For a finite nucleus with a charge density 
distribution q{x) the 60 , 2 , 21 / coefficients in Eq. [^are replaced 
by other coefficients co, 2 , 2 i/ involving integrals of the form 
factor q{x) (see Appendix A). If we approximate the charge 
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density of the nucleus as a uniform sphere the resulting form 
factor takes this form: 


q{x) 


9 sin a; 

—-cos X 

X X 


(4) 


As a further simplification we can approximate q{x) with a 
Gaussian function q = e . The coefficients Cq^ 2 , 2 v can then 
be expressed analitically as simple functions of Z, a and the 
nuclear radius (Appendix A). Unfortunately also this method 
faces numerical limitations that we will discuss in Sect. 0 It 
will anyway be used for small scattering angles as a validation 
for the Monte Carlo approach which we are about to describe 
in Sect. |In] 


III. Monte Carlo simulation of multiple 

SCATTERING 

Multiple Coulomb scattering can be simulated numerically 
using different approaches: 

• detailed simulations: each track is simulated at micro¬ 
scopic level as a succession of connected straight seg¬ 
ments (free-flights), Changes of direction are determined 
by sampling the polar deflection cos 6* from the single 
scattering distribution, and the azimuthal scattering angle 
uniformly in [0, 27r]. The length t of the free-flight is sam¬ 
pled from the exponential distribution p{t) = 

The simulated angular and spatial displacements are 
“exact”, i.e. they coincide with those obtained from a 
rigorous solution of the problem (i.e. by mean of the 
transport equation introduced in the theories of Goudsmit- 
Saudersons El, ii and Lewis 0). Nevertheless detailed 
simulations are in practice not possible computationally 
except for specific cases involving very low energies and 
thicknesses. 

• condensed simulations: the global effects of the col¬ 
lisions during a macroscopic step are simulated us¬ 
ing approximations. Usually the theories of Goudsmit- 
Saundersons El, a, Moliere Q, a and Lewis a are 
used to predict the angular and spatial distributions after 
each step. At energies where condensed algorithms are 
required, the great majority of elastic collisions are “soft” 
collisions with very small deflections. Especially the spa¬ 
tial distributions might present significant dependencies 
on the choice of the macroscopic step which has to be 
carefully studied. 

• mixed algorithms: Hard collisions, with scattering angle 
9 larger than a given value 9s, are individually simulated 
and soft collisions (with 9 < 9s) are described by 
means of a multiple scattering approach. It is clear that, 
by selecting a conveniently large value for the cutoff 
angle 9s, the number of hard collisions per electron 
track can be made small enough to allow their detailed 
simulation. The value of the mean free path between the 
simulation of hard collisions, X{h), is given by = 

a{9) sm9d9. As the fluctuations in the spatial 
displacement after a certain path length are mainly due 
to hard collisions, a mixed procedure (with a suitably 
small value of 9s) yields spatial and angular distributions 



Fig. 1. Distribution of the step length for the GEANT4 mixed simulation for 
2 GeV/c muons over 1 mm of Lead. 


that are more accurate than those from a condensed 
simulation. 

Since we are interested in rare high-g events a mixed 
algorithm approach has been preferred. Furthermore in this 
scheme the inclusion of nuclear effects can be easily taken 
into account when sampling the angular distribution at single¬ 
scattering level. 


IV. The Geant4 Monte Carlo implementation 


The version of Geant4 used for this study is 4.9.6.P02 
with the standardSS physics list option to enable the treat¬ 
ment of single scattering above a certain angular thresholcf] 
The simulation of single-scattering is dealt with by the routine 
G4WentzelOKandVIxSection . cc. The distribution of 
the step length used by the GEANT4 simulation is shown 
in Fig. [T]for 2 GeV/c muons on 1 mm thick lead target. The 
slope of the exponential distribution corresponds to a value of 
^(h) = 0.14 pm such that the single-scattering code is invoked 
on average about 7000 times for each muon crossing the 1mm 
thick lead target. 

The event sampling is done according to the model by 
Wentzel El in the formulation of Qo). The model adopts a 
modihed Coulomb potential to take into account the effect of 
the screening of the electron cloud: V{r) = (Ze^ jr)e~'^^^. 
A slightly modified version of the Moliere screening angle 
x'a = ( 2 ^)^ whose Z and energy dependence is 

parametrised as (cfr. Ha, eq. 23): 


Xa = 


9 2 

a 


ZV3 


1.13-f 3.76 


aZ 

T 


T + Z^f^ 


(5) 


ao being the Bohr radius and r = Ekm/rn. 

If we define p = i(l ~ cos0) the differential cross section 
da/dp in presence of screening is proportional to ^ P~^(Xa + 
p)“^. The inverse transform method leads to the following 
sampling formula: 


(Xg +MS)C(1 - IJ-S) _ 

(Xa + 1) - C(1 - Ps) ~ 1 - C(1 - Ms) 


(6) 


*“.mac” file setting: /phys/addPhysics standardSS. 
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with fis = 5(1 ~ cos0s)- In the multi-GeV regime the 
screening parameter x'a found to be of 0(10“^^) and hence 
the sampled distribution effectively reduces to a Rutherford 
law dN/dO oc (1 — cos0)“^. The single scattering angular 
distribution generated from Eq|^ is then corrected for effects 
related to the atomic electrons, the nuclear recoil and the 
nuclear charge density by introducing a function grej{9) — 
ga{0) X grecifi) X |F/v(0)p. A random number C' is thrown 
between 0 and 1 and the scattering is discarded if C,' > grej- 
The atomic coiTection reads ga{9) = 1-/3^ sin^ | + + 

cos 6) sin I while the recoil correction follows the term found 
in the Mott scattering formula grec = (l + (1 ~ 

In the standard implementation nuclear charge ef¬ 
fects are implemented through a dipole form factor 
|TAr(<?)P = (1 + ^ with g ~ p0 in MeV. The 

nuclear extension is parametrised as Rjsi = 1.27A° giving 
Rn = 5.37 fm for the nucleus. Dipole form factors 

arise from an exponentially decreasing nuclear charge density 
and are thus not appropriate to describe the charge distribution 
of large nuclei. We have then modified the simulation program 
by introducing an ad hoc form factor for Lead based on the 
Saxon-Woods parametrisation of the charge density which is 
known to produce a much more realistic description; 

Psw{r) = Po{l +e~^ (7) 

The Saxon-Woods form factor has been fitted on precise 
electron scattering data available in the literature HD. We have 
used b = 6.647 fm and a = 2.30/(41n3) fm as in lfT2l . The 
dipole form factor predicts a largely overestimated contribution 
of high-g scatterings due to the assumption of a charge density 
being much more concentrated with respect to what is actually 
observed. It is hence only appropriate for small g. The form 
factor |F^'^(g)p has been then extracted numerically by 
taking the squared Fourier transformation of the charge density 
Psw {’>’)■ The analytical dipole formula for |TV(g)P has been 
then replaced by a look-up table representing the Saxon-Wood 
form-factor of in bins of g (100 bins from g = 0 to g = 788 
MeV/c). Values above 788 MeV have been neglected. We have 
also considered the contribution of scattering on single protons 
(Fp) by composing the two effects according to the relative 
charge Z following the procedure of Butkevic et al. (HD): 

\Fiq)\^ = \F^^{q)\^ + ^(1 - \F^^{q)\^)\Fp{q)\^ ( 8 ) 

Equationj^is the same as Eq. 4 of (lUD) neglecting the atomic 
form factor which is only relevant at very small momentum 
transfer. For the proton form factor Fp we have used an 
exponential charge density (dipole) with p{r) = R~^e 
and Rp = 0.71 fm. 

The Lead form factor \F^^{q)\'^ used in the simulation is 
shown in Fig. (solid curve). The inclusion of scattering off 
protons (dotted curve) increases significanlty the probability of 
large-angle scattering as shown in Fig. bottom, as it can be 
seen by comparing the solid (|T'(g)p) and dotted (|Fjv(g)P) 
curved. In summary the sampling is applied using a g).gj(0) = 
3rec(0) X |F(g)p xga{0). At the interesting energies the effect 
of atomic corrections is essentially absent (g^ > 1 for any 




Fig. 2. Top: charge density for the Lead nucleus according to the Saxon- 
Woods parametrization of Eq. ^ with the parameters described in the text 
(dotted line) and for the proton (continuous line). The relative size is given 
by Z = 82. Bottom: the nuclear form factor obtained by the squared Fourier 
transform of psw(i') (|l^]v('?)P. dotted), the proton form factor obtained by 
the squared Fourier transform of Pp(r) (|^{g)P, dashed). The combined 
form factor (|i^((7)P, continuous) as in Eq.l^ 


angle) and recoil introduces a modest correction due to p/M 
being ~ 0.5 - 7%. 

V. Benchmarking the Monte Carlo with 

THEORETICAL PREDICTIONS 

We have compared the result of the Monte Carlo simulation 
of Sect. HY] for muons of 2, 8 and 14 GeV/c momenta and 
orthogonal incidence on a 1 mm thick Lead foil with the 
predictions of the theory by Meyer discussed in Sect. The 
results are shown in Fig. (2 GeV/c), Fig. (8 GeV/c) and 
Fig.|5] (14 GeV/c). The Monte Carlo distribution is shown by 
the histogram while theoretical predictions are shown by the 
dashed curves. Both the expectation for a point-like (red) or 
diffuse (black) charge distribution are shown. It can be noticed 
that the theoretical predictions show a steep drop above a 
certain value which depends on the muon momentum. This 
is understood in terms of numerical accuracy of the program 


(the sums in Eqns. 15 and 18 involve terms with large powers 
of 0). This behaviour prevents using the Meyer theory for 
a quantitative estimate of LAS, as anticipated. Nevertheless 
we note that, despite the approximations used in deriving the 
theoretical curves (Sect. below the critical maximal angle, 
there is a resonable agreement with the Monte Carlo results. 


VI. Benchmarking the Monte Carlo with 

EXPERIMENTAL DATA 

We have tested the capability of the modified Monte Carlo 
(Sect. to reproduce several data-sets available in the 
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2 GeV/c 



Fig. 3. Distribution of the scattering angle for 2 GeV/c muons impinging 
orthogonally on a 1 mm thick lead target according to the prediction from the 
theory of Meyer (lines) and to the GEANT4 based simulation described in 
the text (histograms). Both the cases of a point-like or extended nucleus are 
shown in red and black respectively. The arrow marks the angle corresponding 
to the OPERA signal region (Sect|Vj^. 


8 GeV/c 



Fig. 4. As for Fig. but for 8 GeV/c muons. 


14GeV/c 



Fig. 5. As for Fig. [^but for 14 GeV/c muons. 


p (GeV/c) 

particle 

material 

t (mm) 

^data 

Ref. 

7.3 


Cu 

14.4 

3.1 X lO"* 

El 

11.7 


Cu 

14.4 

8.7 X 10^ 

El 

2.0 


Pb 

12.6 

2.5 X 10' 

El 

0.512 

e 

Pb 

0.217 

~ (1.25 X 10^'‘)/s 

El 

[1, 15] 


Pb 

2.0 

o(W~) 

da, In] 


TABLET 

Experimental data-sets used for the benchmarking of the 

SIMULATION, ridata INDICATES THE NUMBER OF MEASURED PARTICLES. 


literature. The goal of these experiments was to investigate the 
scattering of muons or electrons at high momentum transfer 
to study the nuclear charge distribution. They are hence 
particularly suited for our purpose despite being rather old. A 
certain degree of extrapolation is present with respect to the 
CNGS case due to the use of different energy and particles and 
targets. This is mitigated by the fact that a significant variation 
of these parameters is sampled and that overall the conditions 
are not too dissimilar to the case of interest. A summary of the 
basic parameters of the used data-set is given in Tab. The 
Copper targets are about a factor 3-4 thicker than the OPERA 
case while the momentum is in the OPERA window. The Lead 
data are represented by a target being thinner by about a factor 
ten (low-momentum electrons) and one being about a factor 
five thicker but in the OPERA momentum window. It must be 
noted that however the characteristic quantity of the process 
is the transverse momentum transfer q rather than the incident 
momentum. The ability of the simulation in describing the 
data in this space is a good indicator of its reliability in the 
OPERA region of interest. 

A. Copper data with muons (Akimenko et al.) 

The experiment 111 was performed in 1986 at the IHEP 
accelerator using the HYPERON setup to measure the scatter¬ 
ing of muons of energies of 7.3 and 11.7 GeV/c off a 14.4 mm 
thick (one Xq) Copper target. The incoming muon direction 
was determined with four proportional chambers allowing a 
spatial resolution of 0.7 mm and an angular accuracy of 0.3 
mrad. A magnet was also employed allowing the incoming 
momentum to be estimated with a 3% relative error. Down¬ 
stream of the target, the positions and angular errors were 
reduced by about a factor 2 w.r.t to the upstream section. Lull 
efficiency is claimed up to 50 mrad (qmax = 1.9 fm“^ for 
7.3 GeV/c and 3.0 fm“^ for 11.7 GeV/c). A 1.5 m thick iron 
absorber was used to stop pions. A selection based on a pair of 
Cherenkov counters upstream of the target allowed to achieve 
a pion contamination below 2%. In addition, the requirement 
of single track scattering was applied to reject events with 
nuclear interactions of the pion contamination. 

The scattering in the counters was deconvoluted by using 
data-sets obtained after the target had been removed. The com¬ 
parison with the simulation is performed for 0 > 5 mrad since 
in this region the background contribution can be neglected. 
Real data and Monte Carlo events are normalised to each 
other above this threshold. The results are shown in Eig. 
and Eig. |7] for 11.7 GeV/c and 7.3 GeV/c muons respectively. 
In this case the same procedure used for Lead was employed 
except for the fact that the Saxon-Wood parametrisation of the 
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Fig. 6. Comparison of the Monte Carlo simulation to the 11.7 GeV/c data-set. 
Data have been obtained by digitisation from the original paner lMI . 

Copper nucleus charge density was used. Experimental data 
are shown by the bullets with error bands. The Monte Carlo 
prediction is shown with three options: 1) a point-like nucleus 
prediction (|E"(q)p = 1) (gray). 2) the dipole nuclear form 
factor (red) 3) the improved form-factor of Eq. (yellow). 
In both cases the mixed algorithm described in 11^ has been 
used. The band width in Monte Carlo histograms represent 
the statistical uncertainty. This meaning of symbols is also 
preserved in the following. The disagreement with the point¬ 
like simulation indicates that the nuclear density is in this case 
actually probed. Eurthermore a better description is obtained 
by using the improved form-factor with respect to the dipole 
one. 


B. Lead data with muons (Masek et al.) 

The experiment ifTsl was performed in 1961 profiting of 
the Bevatron accelerator at the Lawrence radiation labora¬ 
tory employing a muon beam with a median momentum of 
(2.00 ± 0.03) GeV/c and a 3.5% spread. The total number 
of muons incident on the apparatus was 2.5 x 10^. The 
lead target had a thickness of 14.4 g/cm^ corresponding to 
1.268 cm (one half inch). Scattered particles were observed 
up to 12° (pt — 400 MeV/c). The muon beam is obtained 
from a (3.5 ± 0.3) GeV/c pion beam with magnetic selection 
(see Eig. 3 of the original paper). 2 GeV/c muons are those 
decaying backward in the tt rest frame. This large difference 
in momentum allows a clean 7r/p separation. The pion con¬ 
tamination on the target is estimated to be of about 3%. By 
using a 107 cm (42”) thick iron absorber downstream of the 
target and an upstream Cherenkov counter the effective pion 
contamination was reduced to 4.9 x 10“®. Emulsion detectors 
downstream of the target were also employed to control the 
pion contamination. 



0 (mrad) 


Fig. 7. Compaiison of the Monte Carlo simulation to the 7.3 GeV/c data-set. 
Data have been obtained by digitisation from the original naper lHl . 


The incoming and outgoing muon directions were deter¬ 
mined by using four identical counter hodoscopes made of 
scintillator bars at positions Aa, Ab, Ac and Ajj (Eig. 
3 of HSl). The acceptance of such an arrangement was 
between 2 and 14°. Each station was composed of 20 vertical 
scintillators each having dimensions (0.95, 2.54, 15.24) cm, 
the smallest dimension being the one in the horizontal plane. 
The downstream stations centers were shifted away from the 
beam direction by 26.2 cm. The distance along the beam 
between the stations was 2.794 m (AB) and 2.007 m (CD). 

The beam spread in the horizontal and vertical directions 
was determined triggering with a small scintillator (S) located 
along the beam axis in between stations Ac and Ajj, rotating 
the A A and Ab stations by 90°. The coincidence rate between 
strips in scintillators in A a and is shown in Eig. 6 of 
the original paper. This bidimensional rate map was used to 
generate the incoming direction of muons in the simulation 
accordingly. The vertical spread was generated assuming a 
gaussian distribution with a ±0.6° r.m.s. This is supported by 
the paper which states a spread of ±0.8° and ±0.6° in the 
horizontal and vertical views respectively. A gaussian fit of 
the 9x distribution obtained from the information contained 
in the hit strip bidimensional distribution (Eig. 6 of ifTSll l and 
the detector geometry gives a spread of about 0.6° not too far 
from the 0.8° quoted in the paper. 

We have reproduced with a simple simulation both the 
angular spread of the beam and the finite-solid angle effects 
introduced by the strip width. This was performed by binning 
the true positions according to the detector geometry described 
above. A simple linear trajectory has been assumed in between 
different scintillator stations. 

The measured quantity is the angle $ which is obtained 
considering the layout of hit strips in the scintillator counters 
which are only measuring the horizontal projection of the 
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|G4S.W. FF 
I G4 dipole FF 
|G4 FF= 1 
Masek et al. 



Fig. 8. Comparison of the Monte Carlo simulation to the Masek et al. data-set. 



scattering angle (single-view readout). The distribution of $ is 
shown in real data (bullets) and Monte Carlo (histograms) in 
Fig.J^ The meaning of the symbols is the same as in Figs. 
andj^The Monte Carlo normalisation corresponds to the same 
number of incoming muons as in the data (Tab. |I]|. In this case 
both the point-like and the dipole form-factor simulations fail 
to describe the measurement while the modified simulation 
is capable of reproducing the experimental points at large 
scattering angles. The marked disagreement in the first bin is 
supposedly caused by the simplified simulation which might 
lack some of the geometrical details of the real apparatus. 

C. Lead data with electrons (Frois et al.) 

The experiment im was performed in 1977 at Saclay 
directing an electron beam with an energy of 502 MeV and 
a 20 /rA current (^ 1.25 x 10^"^ e/s) on a water cooled 
217 ± 2 mg/cm^ thick lead target. Thanks to the large event 
sample the data allowed a precision on the charge density 
at the lowest values of r ever reached at those times or 
equivalently they allowed to determine the nuclear form factor 
at momentum transfers q up to 3.5 fm“^ (cross sections 
down to 10“^° mb/sr). The measured distribution of the 
momentum transfer q spans over 12 decades (Fig. bullets). 
The scattering angles were determined with an accuracy of 
0.05°. The experimental points are compared in Fig. with 
our simulation. In this case an event-by-event weighting was 
applied to the point-like distribution following the form-factor 
of Eq. lousing the integral momentum transverse in the target 
as reweighing variable. This was dictated by the difficulty 
in adopting the single-scattering sampling approach within a 
reasonable amount of CPU time. It was verified anyway that 
the sampling approach correctly describes the experimental 
distribution up to g ~ 1 fm“^. In the high-g region the 
two approaches are expected to be very similar due to the 
prominence of single scatterings. The sharpness of the dips in 


Fig. 9. Comparison of the Monte Carlo simulation to Frois et al. data-set. 


the simulation is an artifact of the event-by-event instead of 
scattering-by-scattering form-factor reweighting. 


VII. Large-angle muon scattering in the OPERA 

EXPERIMENT 


The large-angle scattering of muons from CC interac¬ 
tions represents a background for the t ^ fjL channel in 
the OPERA experiment at the CNGS neutrino beam. The 
kinematic selection described in the experiment proposal na, 
ca defines the signal region for the r —/r search by 
requiring: 

. > 0.25 GeV/c 

. 0^ >2G mrad 

• > 1 and p^ < 15 GeV/c 


Eurthermore the kink is asked to lie within the first two lead- 
emulsion layers downstream of the primary interaction vertex. 
One Vt candidate in this decay channel was observed Cl 
with p^ = (0.69 ± 0.05) GeV/c, 0'" = (245 ± 5) mrad and 
p^ = (2.8 ± 0.2) GeV/c. In the experiment proposal ifTa . 
C3 some data-driven estimates of this process are reported. 
They are summarised in Eig. The red arrow represents an 
upper limit from the non-observation of high-angle scattering 
in 1.5 X 10® mm of muon tracks in the nuclear emulsions 
of the CHORUS experiment. This result was translated to 
2 mm of Lead with the CNGS spectrum yielding oa a 
90% CL limit at 4.1 x 10“®. The blue line corresponds 
to the already described measurement on a Copper target 
(Sect. VI-AI. A limit of 2.3 x 10“® at 90% CL was obtained by 
extrapolating the observation to the case of scattering on Lead 
assuming a scaling going as Z'^pjAp^. The green hatched 
band was obtained from the observation of a single high- 
scattering event in the preliminary phase of a dedicated test- 


beam experiment C3 giving 0.6to;6 x 10 °. In the proposal 
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Fig. 10. A compilation of experimental results on Large Angle muon 
Scattering. Bands represent Icr confidence limits. The plot is discussed in 
the text. 



Fig. 11. Spectrum of CNGS events after selection used for reweighting 
the uniformly distributed momentum distribution used as an input to the 
simulation. 



crossed cells 

with 9 > 20 mrad 

and Pt > 0.25 GeV/c 

n. 

2.22952 xlO^ 

1082520 

115 

weighted 

3.05367 xlO^ 

365195 

179.155 

fraction 

1 

2.4 x \0-‘^lvf^ 

(1.2 ±0.1) x lO-'^/uf^ 


TABLE 11 

Simulation results for the OPERA LAS background estimate 

INCLUDING MOMENTUM RECONSTRUCTION UNCERTAINTIES. 


and later OPERA analyses the background has been finally 
assumed to be 1.0 x with a 50% relative error. 

The contribution of LAS in this kinematic region was 
estimated using the simulation described in |IV| A total number 
of about 1.1 billion incident have been generated with a 
flat momentum distribution between 1 and 15 GeV/c with or¬ 
thogonal incidenc^on the lead-film double cell. The obtained 
distributions have then been reweighted in the incident muon 
momentum in order to reproduce the real spectrum of CNGS 
events. The used reweighting function HD (Fig. ED 
reflects the expectation after applying the full selection chain 
used in the analysis ll20l . 

It must be noted that in the real case of the OPERA exper¬ 
iment the momentum is determined with a certain resolution 
and that the selection cut pT > 250 MeV/c has to be taken 
on the reconstructed variable. The momentum resolution of 
muon tracks depends on whether the muon stops in the target 
or it exits from it. Furthermore if the track does not cross 
the spectrometer the only available determination comes from 
the analysis of multiple Coulomb scattering in the emulsion 
detectors. A very conservative estimate of the achievable 
resolution has been taken by considering a 15% Gaussian 
smearing in the p~^ variable (See f.e. Fig. 91 of mi 

The scattering angle 0^ is defined as the difference in slope 
at the entrance and exit of the lead foil. We show in Fig. [T^ 
the distributions of 0^ (top) and pi^ (bottom) before cuts 
(black), after the angular cut (red) and after the angular and 
Pt cut (magenta). The numbers are break down as the cuts 
are applied in Tab. [I^ The probability of having a LAS event 
in the signal region over two mm of Lead is then found to 

^The angular distribution of muons is such that the average path in 
lead in reality is larger than what assumed here by a few %. 


be: fj^Xs ~ (1-2 ^ 0.1(stat.)) x (yellow band 

in Fig. [T0| ). The effect of momentum smearing convoluted 
with the observed distributions has a sensible effect since it 
increases the background due to migration of events in the 
high-pj’ region by bringing it from to (6.8 ±0.8) x 
to (1.2 ±0.1) X 10-^/7/^<^. 

A rough estimate of the OPERA LAS background can be 
also obtained by considering the data from Masek et al. lfPSll 
discussed in Sect. IVI-BI In that case the muon momentum is 
2 GeV/c which is included in the momentum range considered 
for the OPERA signal region ([1, 15] GeV/c). The thickness 
is 12.6 mm to be compared with the 2 mm of OPERA. The 
projected angle $ is considered in that case. Muons with 
Pt > 250 MeV/c will have approximately $ > 0.088. On 
a sample of 2.5 x 10^ muons from the published binned 
distribution we can estimate the observation of (7 ± 3) 
events. Scaling linearly for the different thickness one can then 
estimate: GeV/c) = (0.44 ±0.20) x 10“^/z/^^ 

which lies even below the Monte Carlo result averaged in the 
[1, 15] GeV/c window (dark red band in Fig. [TOl l. 

A. Contribution of the scattering in other materials 

The decay of the t lepton is also considered when occurring 
in the emulsion films. For this reason we have also evaluated 
the LAS probability in these materials. Due to the high number 
of atomic components the implementation of Saxon-Woods 
nuclear form factors was not attempted in this case. For 
this reason we can take the results as an upper limit of the 
background. It must be also noted that for lighter nuclei the 
use of dipole form-factors, while remaining necessary, is less 
compelling than for the case of Lead. Each film is composed 







































































Fig. 12. Distributions of (top) and 6^ (bottom) before cuts (black), after 
the angular cut (red) and after the angular and px cuts (magenta). 


by a TAC (TriAcetylCellulose, C 28 H 3 gOi 9 ) 200 /iin thick 
plastic base and by a pair of 50 thick emulsion layers 
on both sides of the base. The TAC has been modeled as a 
mixture of Carbon, Hydrogen and Oxygen in the proportion of 
28:38:19. The detailed element budget for the emulsion layers 
ETl in terms of mass fraction is Ag 38.34%, Br 27.86%, C 
13.00%, O 12.43%, N 4.81%, H 2.40%, I 0.81%, S 0.09%, Si 
0.08%, Na 0.08%, K 0.05%, Sr 0.02%, Ba 0.01%. The results 
are: flTs < 1-2 x 10" and /Jff < 0.7 x 10 -^ 1 ^^^. 

VIII. Muon photo-nuclear interactions in the 
OPERA EXPERIMENT 

For completeness we also try to estimate a possible contri¬ 
bution to the T —>■ background allowing for the possibility 
that it might arise from muon photo-nuclear interaction with an 
undetected hadronic remnani^ At high enough muon energies 
{E > 10 GeV), and at relatively high energy transfers e.g. 
when the energy lost by the muon, v — E' — E, is more 
than a few % of its initial energy, the contribution from the 

^We have not considered in detail the possible contribuition from muon hard 
bremsstrahlung since the high energy gamma associated to the large scattering 
of the muon can be vetoed very efficiently in the emulsion detectors (22) 


inelastic interaction of muons with nuclei ||23]| starts to play a 
non negligible role. 

The differential cross section a{E^v) (cm^ g”^ GeV“^) 
can be factorised as a{E,iy) = t/) ll24ll . The first 

function, tp, reads: 

T'(z/) = —Na^^cTjn (9) 

TTl/ A 


with a^N = (49.2 -f 11.1 In -f x 10 cm^ with v in 
GeV The nuclear shadowing effect 1251 is parametrised 
as Aeff = 0.22A-f 0.78AO ®9 following EH (0.65 for Lead). 
The adimensional function $ depends on the fraction v = E/v 
and E as follows: 


^{E, v) = V — l-\- 


1-V+- 


1 + 


2/i^ 


In 


f 1 _|_ 


I Ev (i JE A- 

A A/ 


( 10 ) 


M being the nucleon mass, the muon mass and = 0.4 
GeV2. 

An energy-dependent scattering probability per incoming 
muon, P^m{E), is obtained after integration on v and multipli¬ 
cation by pt = 1.135 g cm“^: P^pj{E) = pt v)dv 

with Vmin = 0.2 GeV and Vmax = E—M/2 ll28l . The result is 
shown by the black line curve in Fig.[^ In the same figure we 
show with the red line the cross section when one restricts to 
the phase space which is relevant for the OPERA background: 
the energjj^ of the outgoing muon (E' = E — v) should be 
1 < i7' < 15 GeV, the scattering angle (sin0 = 
should be 0 > 0.02 rad and E' sin 9 > 0.25 GeV. The result is 
then weighted with the probability density function of muons 
from of CNGS muon neutrinos (Fig. Inland dashed curve 
on Fig. 1^ giving = 2.7 x 10“^ events/(also shown 
as a reodotted horizontal line in Fig. [O] ). It must be noted 
that this has to be considered as an upper limit to the real 
background level since a complete inefficiency in the detection 
of the hadronic system is assumed here for simplicity. Here 


approximate momentum to energy here. 
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we refrain from trying an estimate of the minimal hadronic 
activity which could be detected in the OPERA emulsion 
detectors since the contribution is already small enough for 
practical purposes. 


IX. Conclusions 

The probability for the occurrence of large-angle scattering 
for multi-GeV muons in 1 mm thick lead targets has been 
addressed using a mixed-approach simulation program based 
on the GEANT4 libraries. A realistic parametrisation of the 
Lead nuclear density has been implemented (Saxon-Woods) 
and scattering off single protons has been considered. The 
developed algorithm has been validated by means of theoret¬ 
ical considerations and supported by experimental data from 
the literature. The simulation setup predicts a background for 
the OPERA > Vr analysis in the t —>■ fj, signal region 
of Ilas — 0-f) ^ which is well below 

the values which have been considered so far. Einally we also 
calculated, in the CNGS context, an upper limit on background 
contributions from the muon photo-nuclear process which 
might in principle also produce a large-angle muon scattering 
signature in the detector; < 2.7 x 10“^ events/. 
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Appendix A 

Moliere theory with nuclear eeeects (Meyer) 

Given a target with thickness t, density p, atomic weight 
A, charge Ze and an incoming particle of momentum p with 
charge ze we can define a “characteristic angle”, Xc, as 


_ 

pfic 


AnNAtp 
A ’ 


( 11 ) 


Na being Avogadro’s number. The “screening angle” (xa) 
is defined from the De Broglie wave length of the incoming 
particle (A) and the Thomas-Eermi radius of the atom (a = 
0.53 X 10-1° m/Xi/3) as xl = + S.TOa^) with 

a = Then the B parameter of Eq.j^is determined by 

the transcendental equation: 


B = —\n 


xl eB 


( 12 ) 


where e and 7 here are the Euler and Euler-Mascheroni 
numerical constants. The angle x appearing in Eqns. is 
Y = — ^- 7 = . In the specific case of a 1mm thick Pb foil hit 

XoVB 

by a 10 GeV muon we find Xc = 0.24 mrad, Xa = 1-61 /rrad 
and B = 14.7. The coefficients Ci for a nucleus with nuclear 
charge density q{x) have the following expressions; 


noo 

Co = 0.423-^2/ {1-q{x)){l-e-^^)x-^dx (13) 

Jo 


^http://arohatgi.infoAVebPlotDigitizer 


POO 

Cl = —1.423 — 2 / (1 — q{x))e~^ x~^dx 

Jo 


(14) 


C 21 / = 




112 


q{x)x 


2v-0-x^ 


dx 


(15) 


If a uniform charged sphere is assumed as in Eq. and ap¬ 
proximating the resulting form factor with a gaussian function 
q = 6 “““ , the expression of the coefficients assumes this 
analytical form: 

1 ^ 


Co = 0.423- 1 - 




^ In ^ - I 1 


In 1 + 


Cl = -1.423- 1 - - In l-f ^ 


C2v = 




7 ;+ 


^ —2 


l-l/ 
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